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Abstract 

We consider the simplest bosonic-type S-matrix which is usually regarded as un- 
physical due to the complex values of the finite volume ground state energy. While 
a standard quantum field theory interpretation of such a scattering theory is pre- 
cluded, we argue that the physical situation described by this S-matrix is of a 
massive Ising model perturbed by a particular set of irrelevant operators. The 
presence of these operators drastically affects the stability of the original vacuum 
of the massive Ising model and its ultraviolet properties. 
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1 Introduction 



In recent years, remarkable progress has been witnessed in the study of the scaling region 
nearby the fixed points of the Renormalization Group (RG) of 2-D statistical models. 
Such progress has been possible thanks to new ideas related to the S-matrix theory [I], Q 
and there is now a well-defined program to analyse off-critical statistical models. In 
short, this usually consists in implementing the following steps: (a) a deformation of 
a Conformal Field Theory (CFT) describing the critical point by means of a relevant 
integrable operator which drives the model away from criticality; (b) a consequent de- 
termination of the exact spectrum of the massive or massless excitations and the elastic 
S'-matrix scattering amplitudes thereof; (c) and finally, a reconstruction of the off-critical 
correlators by means of the Form Factor Approach also accompanied by an analysis of 
the finite size behaviour of the theory by means of the Thermodynamics Bethe Ansatz 
(TBA). The scattering theory of the Ising model in a magnetic field together with 
the calculation of its correlators provide an explicit example of a statistical model 
solved away from criticality along the line of the above program. Solutions for many 
other off-critical models have also been found (see for instance [3-19]). Hereinafter the 
statistical models solved according to an S'-matrix program will be simply referred as 
bootstrap models. 

All known boostrap models greatly differ from each other by the nature of their spec- 
trum, for the total number of particles and for the detailed structure of their scattering 
amplitudes. However, looking at them as a whole, they share an important feature, which 
consists in the fermionic nature of their scattering amplitudes. Namely, the S-matrix in- 
volving two identical particles takes the value —1 at the threshold of their s-channel, i.e. 
S aa (0) = — 1 (for the notation see below). This condition was proved to be quite im- 
portant in the further study of their off-shell behaviour. In the Thermodynamics Bethe 
Ansatz calculations, for instance, fermionic-type S-matrices lead to well-defined integral 
equations which - in all known models - provide real solutions for the ground state en- 
ergy, in perfect agreement with the CFT prediction fl4] , [L5[] . On the other hand, in the 
Form Factor Approach, fermionic S-matrices give rise to spectral series of the correlation 
functions with fast convergent properties (see for instance, || |9], [10], [IT], 0). 

Among the set of known solved models the peculiar absence of examples relative to 
an S-matrix of bosonic-type, i.e. an S-matrix which satisfies the condition S aa (0) = 1, 
is indeed intriguing. What is the reason? While a widespread suspicion is that such 
S-matrices are unphysicalQ, a full understanding of their features and also of their inter- 

1 It is worth mentioning that in an example based on the algebraic Bethe Ansatz p3|| , which is a 
different approach to that of this paper, it was shown that the physical vacuum of the theory could only 
be constructed if the particles were of fermionic type. 
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pretation are still missing. 

Deeply related to the absence of models with a bosonic-type S-matrix, there is also 
the question about the role played in the scattering theory by the so-called CDD factors. 
It is well known that they give rise to ambiguity in the S'-matrix of a bootstrap system 
and it seems therefore quite striking that all consistent boostrap models have managed 
to successfully resolve this ambiguity. 

The aim of this paper is to address some of the above questions. In particular, we 
will draw some conclusions about the interpretation of the physical aspects which occur 
in an integrable theory based on a bosonic-type S'-matrix. We will consider the simplest 
non-trivial scattering theory of this kind, consisting of one self-conjugate particle A of 
mass m, with the two-particle S-matrix given by 

SUB) = - Si "^T' Si " Ta - -W) ■ (1.1) 
smh p + i sin ixa 

Unless explicitly stated, a is a real parameter which, for the invariance of the S-matrix 
under a — ► 1 — a, can be taken in the interval [0, 1/2]. In eq. (|1.1|) (3 = fii — p\, where 
(3i is the rapidity variable of each particle entering its relativistic dispersion relation 
Ei = m cosh Pi, pi = msinh/5j. Since the Mandelstam variable s is given by s = (pi + 
P2) 2 = 2m 2 (1 + cosh/?), the threshold in this channel is reached when (3 = and for the 
S-matrix (|1.1[) we have 

n , , f 1 ifa^O ; , , 

Saa(0) = (1.2) 
1—1 it a = , 

i.e. a bosonic-type S-matrix except when a = 0. Notice that the asymptotical behaviour 
of this S-matrix is given by lim^oo S(/3) = —1. As it will become clear in the next 
section, the S-matrix (|1.1|) can also be equivalently regarded as a pure CDD factor 
added to an initial fermionic-type S-matrix. 

A word of caution. In the following we will often use the terminology "QFT" to 
denote briefly and concisely a hypotetical theory underlying the scattering processes. 
The terminology does not automatically mean that this underlying theory is a consistent 
Quantum Field Theory, the peculiar features of which are precisely the object of analysis. 

The paper is organised as follows. In Section 2 we briefly discuss the CDD factors 
and their relation with a perturbed CFT. In Section 3 we discuss the difficulties which 
the bosonic-type S-matrix (|1.1|) pose both in the TBA analysis and in the Form Factor 
Approach - difficulties which prevents its interpretation in terms of a standard quantum 
field theory. Section 4 is devoted to an interpretation of the bosonic-type S-matrix ( |1 . 1|) 
as the one coming from the critical Ising model perturbed by a combination of relevant 
and irrelevant operators. We will argue that this combination induces an instability in 
the ground state of the theory, in agreement with the results of the TBA analysis for the 
bosonic-type S-matrix Finally our conclusions are summarised in Section 5. 
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2 CDD Factors and Deformations of CFT 



To analyse the scaling region around a critical point described by a CFT, one usually 
considers the deformation of the conformal action Ac ft by means of one of its relevant 
operators 

A = A C ft + \ I <&{x)d 2 x . (2.1) 



The resulting action - which describes a RG flow from the original CFT to another fixed 
point - gives rise in the Minkowski space to the scattering processes among the massless 
or massive excitations present away from the critical point |§ . Important simplifications 
occur when the action fl2.1p defines an integrable quantum field theory: in this case, in 



fact, all the scattering processes are purely elastic and can be expressed in terms of the 
two-particle scattering amplitudes S a b(/3). On a general basis, these amplitudes satisfy 
the unitarity and crossing symmetry conditions^ 

S ab ((3)S ab (-/3) = 1 ; 

Sab(P) = S ab (ilT - (3) . 

Possible bound states which occur in the scattering processes are promoted to asymptotic 
states so that the amplitudes satisfy the additional bootstrap equations 

SM = S hd {f5 - iu a bc )S ad {(3 + iul) , (2.3) 

where u = tc — u and iu c ab is the location of the simple pole in S a b{fi) which identifies the 
bound state A c in the scattering channel A a A b — > A a A b . 

The solution of the above equations presents the well-known ambiguity relative to 
the CDD factors. Let S a b(/3) be a set of functions with a minimum number of poles which 
solve both eqs. fl2.2|) and (|2.3j ). Another solution of the first equations (|2~2| ) can easily 
be obtained by multiplying each amplitude S a b(fl) by the so-called CDD factors, i.e. an 
arbitrary product of functions f a {(3) defined in eq. S a b([3) — ► C a b(/3)S a b(P), where 

C a b(/3) = Yla^ faAP)- I n i&ct, each f a (P) automatically satisfies eqs. (|2.2|) and hence 
their product. Additional constraints on such CDD factors may come from the bootstrap 
equations (|2.3| ), which they also have to satisfy. Since the CDD factors do not introduce 
extra poles in the physical sheet < Im (3 < n, the conclusion is that the knowledge of 
the structure of the bound states alone cannot uniquely fix the scattering theory. This 
ambiguity is however not harmless since scattering theories which differ by CDD factors 
have a distinct physical content, in particular a different ultraviolet behaviour^ 

2 For simplicity we consider here theories with a massive non-degenerate spectrum of self-conjugate 
particles. 

3 A well known example is provided by the Affine Toda Field Theories |l7], which share the same 
structure of bound states of some deformed CFT but have different central charges . 
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How can we interpret S-matrices which differ by CDD factors in terms of an under- 
lying action? First of all, the CDD factors do not spoil - by construction - the integrable 
nature of the original S-matrix, simply because they only add extra phase-shifts to the 
original elastic scattering amplitudes. Moreover, they also do not alter the structure of 
the bound states - an infrared property of theory -, because they do not introduce extra 
poles in the S'-matrix. Since they mainly influence the ultraviolet behaviour, it seems 
natural to assume that S'-matrices related to one another by CDD factors are associated 
to integrable actions which, in the vicinity of a fixed point, may differ for insertion of 
irrelevant scalar operators^] rji(x) 

A = A C ft + ^ J §{x)d 2 x + ft J r]i(x)d 2 x . (2.4) 

These irrelevant operators should preserve the set of conserved currents (or part of it) of 
the original action ( |2.1| ) and they should also share the same symmetry properties of the 
relevant operator Q(x). Their presence has however several consequences. In fact, from a 
RG geometric point of view the irrelevant operators t]i(x) are associated to RG trajectories 
which flow into the fixed point, whereas the relevant operators are associated to those 
RG trajectories which depart from it. Hence, the simultaneous presence of relevant and 
irrelevant operators gives rise to a RG trajectory which will pass very closely to the fixed 
point in question but without stopping at it (Figure 1). In its vicinity, its action can 
be parameterised as in ( |2.4| ). However, in contrast with the QFT defined by the action 
(PTT|), the action given by eq. ( |2.4] ) does not automatically ensure the consistency of the 
relative quantum field theory in the entire RG space of the couplings, i.e. on all possible 
distance scales. Obviously, there are some preliminary steps in order to give a meaning 
to the expression ( |2.4|) . 

First of all, the irrelevant operators initially present in ( |2.4| ) will be accompanied by 
an infinite number of counterterms, so that the final form of the action of the theory will 
be actually given by the one of (|2.4| ) together with this infinite tower of higher irrelevant 
operators. This status of art however is not catastrophic as it may appear since the theory 
can still remain predictive in view of the fact that all counterterms can be constrained by 
the request of integr ability EO, EH]]. A real source of problem is another one. In fact, in 



contrast to those theories which are obtained by perturbing a CFT by means of relevant 
operators, where their ultraviolet behaviour is completly under control and ruled in fact 
by the initially CFT, the ultraviolet behaviour of those quantum field theories defined 



as in {2A) will be, instead, a-priori unknown. Since at the short distances the effective 



4 It is worth mentioning that the same line of reasoning has been successfully applied in the context 
of Boundary Quantum Field Theories, with the CDD factors entering the reflection scattering matrix 
put in correspondance with the irrelevant boundary deformations of the theory, see [ fl9| . 
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coupling constants of the irrelevant operators become extremely large, it should not be 
a surprise to find out that one may define a consistent quantum field theory in this 
regime only for special sign of the leading irrelevant couplingFl. This condition is realised, 



for instance, for the roaming trajectory of the Toda Field Theories [^2j which in each 
neighborhood of the fixed points it passes, may be regarded as a RG flow induced by 
a combination of the relevant operator $ = <£> 13 and the irrelevant operator 77 = y^i, 
for an appropriate set of the couplings [24]. On the contrary, another choice of the sign 
of the couplings nearby the critical points may induce an ultraviolet instability of the 
ground state of the theory at certain scale A which will spoil a quantum field theory 
interpretation for distance scales R < A. As we will see, this seems to be the scenario 
associated to the S-matrix ( |1 . 1| ) . 



3 Difficulties with the bosonic— type S— matrix 

A deceptively simple S-matrix affected by a CDD ambiguity is obtained by multiplying 
a scattering amplitude S aa (/3) by —1 (which may be considered as a particular f a term, 
i.e. limA^oo /l+iA(/^))- ^ ^ ne or igi na l scattering theory does not have any bound states, 
i.e. if there are no constraints coming from the non-linear bootstrap equations (|2.3| ). 
there is evidently no obstacle to this multiplication. Consequently, by means of this 
CDD factor we can swap from a fermionic-type S-matrix to a bosonic one and viceversa. 
This is the case, for instance of the S'-matrix of the Sinh-Gordon model, which is given 
by S Sh = /b(/3), where B is a function of the coupling constant of the model |J. The 
result of this multiplication is in fact the bosonic-type S'-matrix (|1 . 1|) . For this reason, 
the model associated to the S-matrix (|1 . 1|) may be referred to as "minus Sinh-Gordon 
model" (mShG). 

There is, however, another way of looking at the S-matrix (|1 . 1|) . It may be seen, in 
fact, as the one obtained by multiplying S = — 1 by a CDD factor f a {(3). S = — 1 is the 
S-matrix of the thermal Ising model |4J and therefore the analysis of this bosonic-type 
S-matrix may be pursued along the considerations presented in the previous section, 
as will be done in Section 4. In the meantime, let us investigate what is the nature 
of the problems which are posed by the bosonic-type S-matrix ( [Lip , regardless of its 
interpretation. We will initially study its finite size and ultraviolet behaviour by means 

5 There is a simple and even trivial analogy of this situation with other quantum field theory problems: 
let us consider for instance the hamiltonian of massive boson H = | ((V$) 2 + m 2 $ 2 ) + A<I> 4 to which we 
add the interaction /i$ 6 . Assume that the solution of the theory at fi = is exactly known. Although 
a perturbation theory can be formally applied for any sign of the new coupling fj,, it is obvious that 
the new hamiltonian will define a consistent theory only for positive values of this coupling since for /z 
negative, no matter how small it may be, there is no longer a stable vacuum. 
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of the TBA approach and we will then proceed to the investigation of its Form Factors. 



3.1 Thermodynamical Properties 



There is a simple way to see that the theory defined by the S'-matrix (1.2) presents some 



problematic aspects from the point of view of a quantum field theory. Let us compute 
in fact the free energy F(mR) of such a theory on a cylinder of length L and radius R 
(with L ^> R) by means of the TBA approach [y]]. As usual, -R" 1 can be interpreted 



as the temperature of the one-dimensional system living on an infinite line, = T. 
As it is well-known, once the free energy F(mR) has been computed, the ground state 
Casimir energy Eo(mR) - parameterised in terms of the effective central charge C(mR) 
as 

E (mR) = -^-C{mR) , (3.1) 
bR 

will be given by 

E (mR) = F(mR) . (3.2) 

In the TBA approach, the free energy F(mR) is initially expressed in terms of the 
functional 

= £(p w ) - isfopM) - M/"(p (r) ) . (3.3) 

The energy term is given by 



= / mcosh(3p {r \p)dp , (3.4) 
the entropy term, in the bosonic case, is given by 

S(p,p^) = Jdp[(p + pV)]og(p + pM)-p]ogp-pW]ogpU] , (3.5) 
and the particle number is expressed as 

AT(pM)= f pU{p)d/3 • (3.6) 



In the above formulas, p^ r \(3) is the density of occupied states (roots) of rapidity (3 
whereas p{(3) is the total density of states (roots and holes), p = p^ + p^ h \ They are 
related each other by the integral equation 

p((3) = £ cosh/3 + J p(P ~ P')p {r) (P')^ , (3.7) 

with the kernel <p(/3) given in our case by 

, n . . d 2 sin 7ia cosh 3 

= lo s s = ■ rr^T ■ 2 — • 

dp sinn p + sin wa 
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Let us assume that the functional (|3.3|) admits a minimum with respect to the distribu- 
tions p(P) and p( r \(3), with eq. ( |3.7| ) acting as a constraint. If this would be the case, the 
partition function Z(L, R) in the presence of a chemical potential p (z = e^ R ) and in the 
limit L — ► oo can be computed in the saddle point approximation and expressed as 



Z(L, R) ~ exp \-mLF{R)} , F(R) = [ cosh/? log(l - ze~ eW ) ^ , (3.8) 

L J J-oo 27r 

where F(R) is the value of the functional Rf(p,p( r ') computed at its minimum and the 
minimum condition is provided by the integral equation satisfied by the pseudo-energy 

/+°° (jff 
<p((3- P')\og(l- ze'^)— . (3.9) 
-oo 2tt 

In terms of e(/3), the distributions p(/3) and p^ r \P) which provide the minimum for the 
functional ( |3.3|) are related each other by 

/ = P^. (3.10) 

The above discussion briefly summarises the basic ideas of the Thermo dynamical 
Bethe Ansatz approach. Let us now analyse how they are actually implemented in the 
case of our bosonic model, starting our analysis from the expression of the free-energy 
obtained in the saddle point approximation of the TBA equations, i.e. from eqs. ( |3.8| ) and 
( |3.9|) . At zero chemical potential z — 1, no matter how small a may be (but never zero), 
the numerical solution of the above equations shows the existence of a critical value of 
the variable R, such that for R < R c (a), the free energy assumes complex values (Figure 
2). It is worth noting the novelty of this situation, which in fact does not occur in all 
other known (fermionic) models, analysed in the past in terms of the TBA equations (see, 
for instance, [[HJ). How can we interpret the occurrence of this critical point in R and 
the complex values assumed by the free energy for R < R c 7 What are the consequences 
of this result? The answer to these questions may be expressed both in physical or in 
purely mathematical terms. 

The simplest physical answer is the following. Since the TBA equations provides the 
full resummation of the virial expansion of Z(L,R), R~ 1 (a) is then identified with the 
radius of convergence of the corresponding series. The existence of this singularity indi- 
cates a non-trivial ultraviolet behaviour of the theory, in particular, it definetly precludes 
its standard interpretation in terms of a Conformal Field Theory deformed by a relevant 
operator^. 

6 The problem with the ultraviolet behaviour of such theory can also be seen without having recourse 
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The singularity at R = R c {a) may be effectively regarded as a sort of Bose-Einstein 
condensation phenomenon. To show this, firstly notice that the partition function ( |3.8| ) 
admits an interpretation as that of a free bosonic gas of exitations \x(P)) Du t with disper- 
sion relations which depend on the temperature fl3fl . In particular, for the energy S(j3) 
of one of these excitations we have £(/3) = e(/3)/R, where e(/3) is solution of the 

integral equation ( |3.9|) . Hence, in a way completely similar to the fermionic case anal- 
ysed in [pql , one can establish the equivalence between the partition function obtained in 
eq. ( |3.8j ) and the one computed in terms of the following series 

oo 

Z(L, J2) = £V (X(A») • • • xWiMPi) ■ ■ ■ x((3n))e-^ )+ - + ^ M) , (3.11) 

n=0 

where the scalar product of the states is computed according to the commutation rules 
of free bosonic particles, with a regularization given by 

[&{P - P')] 2 = ^ cosh(/?) 8{(3 - (3') . (3.12) 



In fact, for the partition function in ( |3.8| ) we have the following expansion in powers of 
(mL) 

Z(L, R) = l- (mL)F(R) + ^£t{F{R)f + ■■■ (-l)-^L( F {R)Y + ■■■ (3.13) 
where F(R) admits the representation 

F(R) = -J2~UR) , (3-14) 

with 



n 

n=0 



I n (R) = [ — cosh/3e- ne(/3) . (3.15) 

J 27T 

On the other hand, computing the partition function by using the expression ( |3.11| ), 

Z(L, R) = l + zZ l + z 2 Z 2 + ■■■ z n Z n + ■■■ (3.16) 
for the first term we have 

Zi = /fW^^/fW^')^^^ (3.17) 
= m L / — cosh(/3) e- £(/3) = (mL) h ■ 



to numerical integration of the above equations. In fact, in consistent theories, the ultraviolet behaviour 
is ruled by the so-called kink configuration of eq. ( |3.9| ) and the corresponding dilogarithmic functions. 
In our case the plateau of the kink configuration, solution of algebraic equation e p i = log(l — e~ £pi ) is in 
fact complex, e p i = ±i^. 



S 



and similarly 



Zi = \{mLfll + \{mL)h ; 

[mLf 4 (mL) 2 (ml) 
^3 - — 3| — ^1 + J i J 2 + —h 



(3.18) 



(mL) 
H ; — J 4 



etc. Hence, it is easy to see that the series ( |3.13|) precisely coincides with that of eq. (|3.16|) , 
the only difference being in the arrangement of the single terms: the terms proportional 
to (mi) in all Z n gives F(R) as their sum, whereas the sum of the terms proportional to 
(mL) n appearing in all Z n gives rise to the higher power (F(R)) n . 

In the light of these considerations, it is then natural to expect that the partition 
function ( |3.8j ) will have a singularity at that value of the temperature i? ( 7 1 (a) where 
the energy of the excitation (at zero rapidity) vanishes (Figure 3). When this happens, 
the contribution coming from the average occupation number with 8 = becomes as 
important as the entire series, producing a singularity in the free energy which resembles 
that of the usual Bose-Einstein condensation. One should not be surprised that such 
effect occurs in this low dimensional system due to the entanglement of the energy 8(0) 
of its excitation with the temperature itself of the system. In conclusion, from the TBA 
analysis one infers that in our bosonic-type model there is a scale of distance A = R c 
below which there is a change in the ground state of the theory. Obviously it is not 
surprising to find out that the evaluation of the series beyond its radius of convergence 
provides complex values, although no direct physical meaning can be assigned to them. 

From a purely mathematical point of view, the occurrence of complex values in F(R) 
simply indicates the failure of the minimization procedure applied to evaluate the free- 
energy ( |3.3|) . Namely, in our case, it is no longer true that the functional ( |3.3| ) admits 
a minimum with respect to the (positive) real distributions p((3) and p^((3) for any 
value of R. As a matter of fact, the real minima of the functional ( |3.3| ) disappear for 
R < R c , becoming complex. Hence the value of the functional (|3~3|) computed at its 
minima, i.e. RF(R), takes for R < R c complex values whereas the functional f(p,p^) 
itself becomes instead unbounded from below for this range of R. This implies that the 
saddle point procedure employed by the TBA is no longer justified^. However, if we still 
insist on enforcing the validity of the relationship (|3.2| ) between the free-energy and the 
Casimir energy of the theory, we should conclude that the divergence of the free-energy 
for R < R c implies the corresponding divergence of the central charge of the theory. As 



7 This circumstance is analysed in full detail in Appendix A for simplified version of a TBA system 
based on a particular bosonic S'-matrix. 
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we will show in the next section, the same conclusion is also reached by using the Form 
Factor approach. 



3.2 Form Factors 

Let us now analyse the computation of the Form Factors associated to the S'-matrix ( |1 . 1| ) . 
As we will see, this calculation presents some new and distinct features with respects the 
analogous calculation for fermionic-type theories. 

The definition of the Form Factors of a local hermitian scalar operator O is given as 
usual by 

F°(Pi, ...,/3 n ) = (0\O\A(f3i), A(p n )) . (3.19) 
These quantities satisfy a set of functional equations || [7| 

F n (Pli ■ ■ ■ i Pit Pi+li • • • 5 Pn) = F n (f3 1 , . . . , fii+i, Pi, . . . , P n )S(Pi — Pi+i) , 

n 

F n (Pi + 2m, . . . , /?„_!, (3 n ) = F n ((3 2 , . . . , /3 n , A) = I] S(J3i - Pi)F n {(3 x , ...,(3 n ) (3.20) 

i=2 

jm0-i3)F^ +2 + i'K,i3,i3 1 ,i3 2 ,...,i3 n ) = i (l - f[S{/3 - A) j F° (A, • • • , Pn) , 

which closely resemble those of the Sinh-Gordon theory |26| . 

For scalar operators, the FF depend solely on the rapidity differences /% = $ — f3j. 
As it is evident from the residue equation in ( p.20| ), the chain of FF with even and odd 
number of particles are decoupled from each other. For simplicity, in the sequel we will 
consider only FF of even operators under the Z 2 symmetry of the model (i.e. those 
which have only non-vanishing FF on an even number of external states). They can be 
correspondingly parameterised as 

FniPi, ■ ■ ■ , Pn) = H n gw( 7'\-;' Xn) TT , (3-21) 

K) 2 y Xi + Xj 

where Xj = e^, H n is a normalization constant and F min (j3) is a function which satisfies 
the equations 

F min (P) = S(P)F min (-P) , (3 22) 

F min (iTc-{3) = F min (iir + P) . 

Q n (xi, . . . , x n ) are symmetric polynomials in the variables xi, . . . , x n of total degree 
t n = f (2n — 1). They can be expressed in terms of the so-called elementary symmetric 
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polynomials given byQ 



a 



(n) 



% i ^ % 2 ' ' ' 



(3.23) 



il<i2<---<ij 

Equivalently, they can be obtained in terms of the generating function 



Y[(x + Xi ) = Y,X n -1<T$ n \ Xi ..r, 



Xn) 



(3.24) 



i=l 



fc=0 



Although at a formal level, the above formulas of the FF appear identical to those of the 
Sinh-Gordon model (except for the extra factor (a n )^ in the denominator), a substantial 
difference between the two models occurs in the properties of F min {{3). As we will see, 
this difference has far-reaching consequences. 

The two-particle minimal FF of our bosonic-type model may be taken as 



F min (/3,a) = exp 
which can also be expressed as 

oo 

F min (P,a) = Af(a)U 

k=0 



dxcosh [§ (1 - 2a)] . 



x cosh | sinh x 



sin 



2^ 



(3.25) 



T [k + \ + % + {k + 1 - ± + ^ 



2t 



rffc + |-f + f r (fc + i + f + g 



(3.26) 



where (3 — vk — (3. With F min (iii, a) = 1, we have 

MM fr/' r(fc + i + S)r( fc + j-f) 
N{a) ll^ r(t + 1 _, )r(i + i + f) 

For a 7^ 0, this function has a finite value at the threshold [3 = given by 

dx cosh [|(1 - 2a)] sinh f 



(3.27) 



F min (0, a) = J 7 (a) = exp 



x 



cosh 2 f 



(3.28) 



which is in stark contrast with the value F min (0) = of fermionic-type theories, which 
is forced by the condition S(0) = — 1 of their S'-matrix. The asymptotical behaviour of 
Fmin{(3,a) is ruled by 



lim F min ({3, a) ~ i \J2 sin7rajF(a) exp 



2 



(3.29) 



8 In the following we will discard the upper index n in the definition of Ofc, since the total number of 
variables involved will be clear from the context. 
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Finally, in the limit a — > 0, F min (j3, a) reduces to 

F min {/3,0) = . (3.30) 

sinn ^ 

In order to implement the recursive equations associated to the last equation in ( |3.20| ) it 
is useful to consider the functional equation! 

2m 2 M 2 (a) 

FminiP + wr, a)F min {(3, a) = — — — , . . (3.31) 

sinn p + 1 sm na 

The residue equation in (|3.20|) implies that the polynomials Q n (xi, ■ ■ ■ ,x n ) in eq. (|3.21|) , 
satisfy the recursive equations 

Qn+2{ X, X\, . . . , x n ) = x D{x, X\ y . . . , X n ) Q n (%l, • • • j 2*n) j (3.32) 

where 

i- fc - p - q sin[7ra(p - g)} jn)jn)„ Mj(n) _ 

(3.33) 



^(X^!, . . . |JBn ) = 2 (-l) 1 ^ V ( _ 1 )^ % 4n-/- fe - P - g ^L^ ~ W UU 

sin7ra ^ 

{Z,fc,p,g}=0 



In writing eq. (|3.32|) we have posed for convenience 



sm7ra 

H n -\-2 = H n - — -— . (3.34) 

(2isin7ray r (a)) n 

Since T>( ) contains four elementary symmetric polynomials and each of them 

is linear in the individual variables Xi, the partial degrees p n of the polynomials Q n 
satisfies the condition p n+2 < 4 + p n . 

After these general considerations, let us now attempt to solve the above recursive 
equations for a particular but significant field. Namely, let us assume that there is an 
operator in the theory which plays the role of the trace Q(x) of the stress-energy tensor 
T IMU (x). If such an operator exists, its Form Factors have some distinguished properties 
which may be useful in their explicit determination. First of all, its normalization will 
be given by 

F e {in) = (A(p)\e(0)\A(j3)) = 2vrm 2 . (3.35) 

Secondly, from the conservation law d^T^x) = 0, the polynomials Q n (xi, . . . ,x n ) can 
be factorised as Q n (xi, ■ ■ ■ ,x n ) = (Ticr n -iP n (xi, . . . ,x n ). Moreover, assuming that this 
operator is even under the Z 2 symmetry of the model, its only non-vanishing FF will 
be those with an even number of external particles. Finally, the two-point correlation 

9 Plugging j3 — in (3.31), one obtains the relation T(a) = 27T J^- — . 
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function of Q(x) enters a sum-rule which permits evaluation of the difference of the 
central charges of the theory by going from the large to the short distance scales BO 

AC = C uv - C ir = ^ J rf 2 x|x| 2 (0|6(x)e(0)|0) = J dp alp) , (3.36) 

where ci(fi) may be directly expressed in terms of the FF of 

*M = ^g^^l^ M\> (3-3T) 

x S(^^m sinh Pi) 5 (^^m cosh f3i — fi) . 

i i 

Since at the large distance scale the theory is massive, we will have Cj r = 0, and the 
above sum-rule should provide a determination of the central charge C uv of the theory 
in its ultraviolet regime. The first approximation of the above sum-rule is obtained by 
the two-particle contribution 

Aci2){a) = ir |Fmm(2/? ' a)|2 • (3 - 38) 

With this information on the FF of Q(x), let us proceed in the computation of its first 
representatives. Since in the limit a — > 0, our bosonic-type 5-matrix reduces to the one 
of thermal Ising model - for which the only non- vanishing FF of Q(x) is the one relative 
to the two-particle state 

F^(P) = -2nmHsmh^ , (3.39) 
we are forced to take as Q 2 {x\, ^2) 

Q 2 {x u x 2 ) =<7i(<7?-402) , (3.40) 

and if 2 — —\rn 2 . These two quantities play the role of initial values for the recursive 
equations ( |3.32| ). 

Notice that by inserting the two-particle FF in ( |3.38| ), the corresponding quantity is 
an increasing function of the parameter a, which takes its minimum value AC*- 2 -* = | 
at a = (as in the thermal Ising model) and reaches its maximum AC^ 2 ) = 0.876... at 
a = I (see Figure 4). This increasing monotonic behaviour of AC' 2 ' (a) is in contrast with 
the decreasing monotonic behaviour presented by the same quantity for fermionic-type 
theories (see for instance the Sinh-Gordon model |26|). 

Using now the factorised form Q 4 (xi, . . . , x 4 ) = (Ti(T^P^{x\ : . . . , x&) and the above 
expression for Qz(xx, 22), the recursive equation satisfied by P^xi, . . . , £4) becomes 

P 4 (-x, x, x u x 2 ) = -4 [x 8 - x 6 al + x 6 cr 2 + x A a\a 2 - x A a 2 2 - a; 2 ^)] (a 2 - 4<t 2 ) . (3.41) 



13 



In the linear space of the symmetric polynomials with total degree t = 10 and par- 
tial degree p < 5, the above recursive equation does not uniquely fix the polynomial 
P^xi, . . . , £4), which in fact admits a three-parameter family of solutions 

P 4 (xi, . . . , x 4 ) = Ai (of 0-304 + 0i0g - afa 2 al) + A 2 {o-\a\o 4 + 0^03 - 010303) + 

+v4 3 (0^04 + 0^0 4 - 01020304) + 32 010 2 0304 + (3.42) 
+ 160204 — 40 x 020 3 — 4010^03 — 640 2 4 z . 

This arbitrariness in the four-particle FF of Q(x) may be partially reduced by imposing an 
additional condition on its FF, namely that they fulfill the "cluster equations" Jl5, |28], |29 



lim F 4 e (ft + A, ft + A, ft, ft) = 77^ if (ft, ft) if (ft, ft) . (3.43) 

By using the cluster properties of the elementary symmetric polynomials 0^, as deter- 
mined in |y|, together with the asymptotical behaviour eq. (|3.29| ) of F min (P,a) and the 



vacuum expectation value (9) = 2 ™ ; the cluster equation ( |3.43| ) provides the following 
conditions on the constants A\ and A 2 

A 1 = -1 , A 2 = -4 . (3.44) 

However, the constant A 3 cannot be fixed in this way, since all terms proportional to it 
in ( |3.42j ) are sub-leading in the above limit Q3.43j ). This arbitrariness of the FF of ©(a;), 



which persists for higher FF, shows that in the mShG model the 5-matrix alone cannot 
uniquely fix the matrix elements of one of its significant operators. This arbitrariness 
in the FF of 0(x) should be contrasted with their unique determination in all known 
examples of fermionic-type S-matrixQ (see, for instance p, [| [T^, f26fl). 



Concerning the c-theorem sum rule, some conclusions can be also drawn despite 
the arbitrariness present in the FF of Q(x). Namely, let us assume that the arbitrary 
constants present at each stage of the recursive equations could be fixed according to 
some additional principle. What will then be the generic behaviour of the series ( p.37[ ) 
entering the c-theorem? For the finiteness of the value assumed by the FF at all particle 



thresholds, an estimate of the integral ( |3.36|) can be provided by only considering the 



contributions at all thresholds. To do so, let us first observe that the sum (|3.37|) - apart 
from the prefactor 12//i 3 -, is an integral of the square of each FF integrated on the 
phase-space f2 2 fc(/u) of 2k particle, defined as 

= I T-^TT " " " 7^%^ - El E ^ ^ + • ■ 'J*) • ( 3 - 45 ) 

J {Z7T)2hi (2iT)2h 2k 



10 It is easy to check that also higher Form Factors of Z 2 odd operators (as the one which creates the 
particle A of the mShG model) cannot uniquely fixed. 
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Near the threshold this quantity may be expressed as 



flafeM - ; —7 ty (A* - 2ibm)"~ 3 . (3.46) 

Concerning the value of the F 2 k(/3i, . . . , fak) a t their threshold = 0, for all z and 
j), we can split it into two terms: the first one is the product of H^u together with all 
(0), given by 

(sin7ra) fe_1 , . ..i^fr-il nm 2 , 9 x k ( T{a) \ k 

H21 J, rwi n (J 7 (a ))* (2fc 1} = (2isin 2 7ra) ^— . 3.47 

(2isinvraJ !r (a)) fc ( fc - 1 ) v v " 2sin7ra v ; \2ismna) v ; 



For the remaining term, we can pose 

Qln\ x li • • • 5 %n) 



°2n Yli<j( X i + X Jy 



fe=0 



^ N fe(2fc-l) 



(3.48) 



An estimate of the ratio p — | may be obtained as follows. From eq. ( |3.24| ), the numerical 
values assumed by the elementary symmetric polynomials o\,- at /3j = (denoted by it,) 

/ 2A; \ 

is given by the binomial coefficient a 3 = ] . The maximum of these values is for 

V i J 

( 2k\ 

cjfc = I ^ 1 ~ A k . At the threshold, in the space of symmetric polynomials of total 

degree k(Ak — 1) and partial degree 2(2k — 1), the term which is expected to dominate 
in the limit k — > 00 is given by [(d'k) Ak ~ z &2k\ — 16 2fc2 ~2 fe , so that £ ~ 16 and the ratio p 
is therefore expected to be always larger than 1. 

Let us now apply the above considerations to the c-theorem sum-rule. We will approx- 
imate the integral ( |3.36| ) by using the "mean theorem" , applied in each interval between 
all successive thresholds. Apart from some constants, the series entering eq. ( |3.36| ) has 
the following behaviour^ 

X 2k \a)Z 2k (a) 



where 



AC^y ± y T rr - (3-49) 
^{2k)\V{k-\)k 1 ; 



X(a) = -^p 2 , (3.50) 
2 sin ixa 



Z(a) 



2 sin 2 71a 



The nature of the series ( |3.49| ) is obviously controlled by the parameter X(a): since the 

2 sin na 



ratio 2 ^inla entering X(a) is always larger than 1.6.. ( see Figure 5) and p is expected to 
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Eq. (j3.49|) holds for a =/= 0. For a = 0, the c-theorem gives AC = h, as in the massive Ising model. 
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be larger than 1, we conclude that, due to the finite value assumed by the FF at the higher 
thresholds, the series Q3.49p is always divergent. This divergence should be regarded as a 
manifestation of the pathological aspect of the quantum field theory associated to mShG 
model in its ultraviolet region, as we already learnt by the TBA approach. 

4 Irrelevant Deformations of the Ising Model 

With the above insights to the nature of the problems presented by the bosonic-type S- 
matrix ( pL.l|) , let us come back to the considerations of Section 2 with the aim to identify 
an underlying action for such an S— matrix. We will take the point of view in which f a {(3) 
in eq. fll.ip is seen as a CDD factor for the fermionic S'-matrix S = — 1. The latter is 
identified as the one of the thermal Ising model. Hence, we have to look for irrelevant 
operators in this model which both preserve its original integrability and its Z 2 symmetry 
such that the extra phase-shift which they induce matches with the CDD factor f a (/3). 

In the conformal Ising model there are few operators we can play with: in fact, we 
have the operators of the conformal family of the identity I (which includes the stress- 
energy tensor T), those of the conformal family of the energy field £ and finally those 
of the magnetization field a. The operators of the even sector of this model, i.e. those 
of identity and energy families, can be written in terms of local expressions of the chiral 
^(z) and anti-chiral ^f(z) components of a Majorana fermionic field. In particular we 
have E(z, z) = i^f (z)^(z) and, for the analytic and anti-analytic part of the stress energy 
tensor, T(z) = \ : ^>d^> : and f = \ : : . 

The original fermionic-type S-matrix S — — 1 originates from the CFT action per- 
turbed by the relevant operator of the energy field 

A = A C FT + m J £(x)d 2 x . (4.1) 

In the euclidean space it may be written as 

A = — / (VOW + f dtf + 2im¥#) d 2 x . (4.2) 

27T J 

Let us analyse which operators we can add to this action in order to have an extra phase- 
shift in the S'-matrix. None of these operators can be of the magnetic field family as they 
would spoil both its Z 2 symmetry and its integrability 0, |32|. They cannot be either 
descendent fields of the energy field, since they are all quadratic in the fermionic field: 
their insertion in the action ( (4.2|) will not induce scattering processes among the fermion 
particles but will only change their dispersion relations. Hence, the first irrelevant field 
which can be introduced in (|4. 1| ) in order to have extra scattering processes is TT . This 
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operator is in fact quartic in the fermion fields and the new action becomes 

A = — I (^d^ + fdf + 2im^) d 2 x + [ dMd* d 2 x + ■ ■ ■ (4.3) 

2n J tt m J 

where g is a dimensionless constant. With the insertion of this operator, the theory is still 



integrable at the lowest order [|33| although A in ( |4.3j ) should be considered at this stage 
as an effective action. The coupling constant g can be related to the parameter a entering 
the S-matrix (|1.1|) as follows. Since the operator TT is expected to be effective at very 
high-energy scales, this suggests matching the lowest order coming from the perturbation 
theory of ( |4.3j ) with the S-matrix ( |1 . 1|) , both computed in their high-energy limit. In 
this kinematical regime the two particles involved in the scattering can be effectively seen 
as left and right movers. Their momenta p and q may be parameterised in terms of a 
new rapidity variable 6 as p = me 01 , q = —me' 02 , so that the Mandelstam variable s is 
expressed in this limit as s = 2m 2 e dl2 . The lowest order in g of the scattering amplitude 



in this regime may be computed as in [20| with the result 



S{s) ^ -1 + 2ig -i- + • - • (4.4) 
m z 

Let us now compare this expression with the one from the expansion of the S'-matrix 
( |1 . 1| ) around the point /?i2 = oo. By setting 1/ sinh/3 12 = e dl2 as a parameterization 
around /3 12 = oo, the S'-matrix ( |1 . 1|) may be written in the vicinity of this point as 

s _L I rn 2 

S(9 12 ) = sn™. _ _ 1 + 2i s i n7 ra — + ••• (4.5) 

sin na 

so that we have the identification 

g ~ sin ira ~ a (4-6) 

i.e. g is a positive quantity. It is important to note that the sign of the above coupling 
constant is opposite of the one relative to the roaming model nearby the Ising fixed point 
||20[1 . This difference in the sign of the coupling relative to the first irrelevant operator 
seems therefore responsible for the quite distinct ultraviolet behaviour of the two theories: 
the roaming model has in fact a smooth ultraviolet behaviour controlled by the nearest 
tri critical Ising point (with central charge C = ^) whereas the QFT associated to our 
bosonic S-matrix develops instead an ultraviolet instability. This situation evidently 
resembles the two different behaviours which originate from the change of sign in the 
example of footnote 5. 

Let us finally add a few comments on the full hamiltonian of the theory. As discussed 
in Section 2, the action (|4.3j ) will also contain higher derivative terms ^d n ^ i &d n ty , so 
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that in general we can write 

A = — I (*<9* + tfdtf + 2zm**) d 2 x + -H—z I tyd^tydty d 2 x + (4.7) 
2ir J % 2 m 2 J 

oo „ 
n=2 J 



The a-dimensional constants a n can be in principle determined by the integrability con- 
dition of the theory, i.e. by matching the expansion of the .S-matrix in power of s/m 2 
with the Feynman graphs which originate from ( |4.8j ). For the presence of the infinite 
higher derivative terms in the action, the corresponding hamiltonian (with respect to the 
Ising fixed point) will be then a non-local one. In the Majorana basisQ denoted as before 
by (ty(x, t), *&(x, t)) , it is convenient to write the Hamiltonian in terms of a dimensionless 
kernel V(z), (z = mz), as 

H = J : [-i^d x ty + Wd x ^ + 2im^] : dx+gm 2 J J : ^ : (x)V(x-y) : : (y)dxdy 

(4.8) 

in such way that the coupling constants of the higher derivative terms are nothing but 
the higher moments of V(z). Therefore, once the coefficients a n were known, the kernel 
V(z) can be recovered in terms of an inverse Mellin transformation. The advantage 
of expressing the hamiltonian in this way is twofold. First of all, it shows that the 
hamiltonian is actually renormalisable albeit obviously non-local. Secondly, it provides 
an easy argument to show that the non-local term is the one responsable for the instability 
of the theory. This argument is based on the simple hypothesis that f V(z)dz > 0. Let 
us consider in fact the Hamiltonian (|4.8|) in the Hartree-Fock approximation and let us 
define an order parameter A = — i(: :). For field configurations which do not vary on 
x, the value of the Hamiltonian (|4.8|) for unit volume is given by 

— = mA - gvA 2 , (4.9) 

where v = J'V(z)dz > 0. When g = 0, we can always subtract from H/L an infinite 
constant (due to the Dirac sea of the fermion) such that we can restrict to the positive 
values of A and the minimum of H/L is therefore obtained for A = 0. This subtraction 
is equivalent to a selection of the vacuum state for the free theory. However, when we 
switch on a positive value of g, the quantity H/L becomes once again unbounded from 
below, so that there would be an instability of the vacuum originally selected by the free 
theory. Hence the above hypothesis permits to understand in an easy way the nature of 
the problem posed by our bosonic S*-matrix, i.e. the instability of the vacuum state of 
the underlying microscopic Hamiltonian. 

12 In this basis 7 = [ | and 7 1 = 

I -i 
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5 Conclusions 



In this paper we have discussed the physical properties of the simplest bosonic-type S- 
matrix with the aim to clarify the nature of its pathological behaviour, in particular at 
its short-distance scales. The TBA approach has shown that the difficulties consist in an 
instability of the vacuum state of the theory which occurs at a certain distance or energy 
scale, i.e. in the failure of the usual saddle point approach to minimise the free-energy. 
This quantity in fact for R < R c becomes in our case unbounded from below. Even in 
the absence of the above method, the peculiar physical behaviour of the QFT associated 
to a bosonic S-matrix may be inferred by the computation of its Form Factors. We have 
seen in fact that even with some stringent constraints imposed on the matrix elements 
they are intrinsically largely undetermined. Moreover, the spectral series which originate 
from them are generically divergent. 

Even after the thorough analysis of this paper, one could be still surprised that an 
innocent —1 in front of an S-matrix is the source of all the pathological features of 
the QFT associated to such an S'-matrix. The fact is that the S'-matrix, even for the 
integrable models where it takes a particularly simple form, is still the final result of an 
infinite resummation of all microscopic processes dictated by an Hamiltonian. Therefore 
while it is in general relatively simple to decide whether or not an Hamiltonian may 
give rise to a consistent Quantum Theory - in terms of stability of the vacuum and the 
relative excitations thereof - (once again, the simple example of footnote 5 is particularly 
instructive), the model analysed in this paper shows on the contrary that it may be 
in general "algorithmically" difficult to trace all this information back and to infer the 
consistency of a theory starting from the end, i.e. from the knowledge alone of the 
S-matrix. 

Several questions arise in relation to the observations of the previous sections. For 
instance, it would be interesting to analyse the QFT associated to more general CDD 
factors and characterise those which have a smooth interpolation from the short to the 
large distance scales, i.e. those ones which identify consistent QFT. Moreover, actions 
which are initially obtained by a perturbation of irrelevant fields contain an infinite series 
of corrections due to the higher dimension operators and hence they are generally non- 
local. In this respect, it would be useful to develop a powerful method for determining 
in an efficient way all the coupling constants of these operators and to have a criterion 
to identify those non-local theories which give rise to integrable models. 
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Appendix A 



Let us consider the calculation of the free-energy of a bosonic type S'-matrix with a 
kernel (p(/3) given by 

<pV3) = 2k6(/3) , (A.l) 

i.e. with a phase-shift defined by 



-ilnS(P) 



-7T if /5 < 

if = 
+ 7T if /3 > 



(A.2) 



In virtue of the locality of the kernel flA.lj ), the total density of levels p(P) is simply 
related to the root density of levels by 



p(P)=a(P)+p^(P) 



where 



in 



a(P) = — cosh P . 



The functional of the free-energy (at zero chemical potential) is given by 



f(p,P {r) )=£(p (r) )-is( P ,p^) 



R 



where 



and 



= / m cos h Pp {r \p) dp 



(A.3) 
(A.4) 

(A.5) 
(A.6) 



S( P ,p {r) ) = J dp [(p + p^)\og(p + p^)-p\ogp-p^ log p^] . (A.7) 

By plugging the expression (|A.3|) for p(P) into eq.( |A.5| ) and minimizing with respect 
to the distribution we have the following condition for the density p^ which 

minimizes the functional ( |A.5|) 

, 2 



(a + 2pW)' 



c 



pM (a + 2pM) 2 
where 

C = e rnRcoshO _ 

For any value of a, the left hand side of eq.( A.8 ) as a function of the positive values of 
is always larger than 4. Hence a positive real solution does not exist for all values of 
R but only for those for which it is verified the condition C > 4. In this case the solution 
is given by 



C-4 



(A.10) 
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Correspondingly 



P{0) 



C 



C-A 



(A.ll) 



Substituing these expressions into (|A.5|) , the value of the functional at its minimum is 
given by 

^(l + Vl-4e- 



-mr cosh f3 



F(R) = m J cosh/? In 
Hence, within the saddle point solution of the TBA equation we have 

E Q (R) = ^G(mR) , 



(A.12) 



where 



G{x) 



dtt 



In 



- (l + Vl - Ae~ xt ^j 



(A.13) 



(A.14) 



This expression becomes complex for R < R c = — . R c is the value where the first 
instability of the functional (|A.5| ) shows up, due to the distribution modes p^(P) and 
p(P) computed at (3 — (for R < R c additional instabilities are induced by other modes). 
To see this, first of all notice that, in view of the algebraic equation ( |A.3| ), we have a 
complete decoupling of the contributions due to different (3 in the free energy. Let us 
consider then the term f° in ( |A.5| ) due to the two distributions computed at (3 = 0. By 



using eq. (|A.3|) , with the notation p 

1 



p( r )(0) and a = ^ we have 



/ (P,R) = P~ ft K a o + 2p) ln(a + 2p) - (a + p) ln(a + p) - plnp] 



(A.15) 



This expression, as a function of p, is plotted in Figure 6 for R > R c , R = R c and 
R < R c respectively. For R > R c f° admits a minimum for a positive real value of p 
and the function is asymptotically positive. Hence for these range of R the saddle-point 
approximation results valid. At R = R c the minimum has moved at infinity although the 
value of the function f° has remained finite. For R < R c the function f° does not have 
any longer a minimum and its values are unbounded from below, causing an instability 
in the corresponding free energy. 
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Figure Captions 



Figure 1 . Renormalization Group flow which passes by a fixed point along the direc- 
tions of the irrelevant and relevant fields 77 and $. 

Figure 2. a . Plots of the central charge c(mR) = —^F(mR) for the mShG model at 
a particular value of a, for the thermal Ising model and for a free bosonic model. 
Below mR c (mR c ~ 1 in the figure), the central charge of the mShG model assumes 
complex values. 

Figure 3 . Profiles of the pseudo-energy e versus (3 for different values of R. R c corre- 
sponds to the value when this function hits the origin. 

Figure 4 . Two-particle contribution to the c-theorem versus the parameter a of the 
^-matrix. 

Figure 5 . Graph of J 7 .^ versus a. 

' - 2 sin wa 

Figure 6. a The free-energy as a function of (0) for R > R c , when there is a minimum. 

Figure 6.b The free-energy as a function of p^(0) for R = R c . The minimum has 
moved at infinity. 

Figure 6.c . The free-energy as a function of /^•'(O) for R < R c when there is no longer 
a minimum. 
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